% impedance rotation
% to verify KKRoy's idea
% latest date - 14.5.3

H = [ -0.5690 - 4.3300i -6.6891 - 4.5159i;-0.4616 - 2.2812i 0.5978 + 0.3821i];

for i = 1:360,
   
   theta(i) = i*pi/180;
   
   Z(i,:,:)=H*[cos(theta(i)) sin(theta(i));-sin(theta(i)) cos(theta(i))];
end;

subplot(221);
polar(theta,abs(Z(:,1,1))');
hold on;
polar(theta,abs(Z(:,2,1))','r-');
polar(theta,abs(Z(:,1,2))','g-');
polar(theta,abs(Z(:,2,2))','k-');
subplot(222);

plot(real(Z(:,1,1)),imag(Z(:,1,1)));
hold on;
plot(real(Z(:,2,1)),imag(Z(:,2,1)),'r-');
plot(real(Z(:,1,2)),imag(Z(:,1,2)),'g-');
plot(real(Z(:,2,2)),imag(Z(:,2,2)),'k-');
